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Abstract 

The role of polymers additives on the turbulent convective flow of a Rayleigh-Taylor system 
is investigated by means of direct numerical simulations (DNS) of Oldroyd-B viscoelastic model. 
The dynamics of polymers elongation follow adiabatically the self-similar evolution of the turbulent 
mixing layer, and shows the appearance of a strong feedback on the flow which originate a cut off 
for polymer elongations. The viscoelastic effects on the mixing properties of the flow are twofold. 
Mixing is appreciably enhanced at large scales (the mixing layer growth-rate is larger than that 
of the purely Newtonian case) and depleted at small scales (thermal plumes are more coherent 
with respect to the Newtonian case). The observed speed up of the thermal plumes, together 
with an increase of the correlations between temperature field and vertical velocity, contributes to 
a significant enhancement of heat transport. Our findings are consistent with a scenario of drag 
reduction between falling and rising plumes induced by polymers, and provide further evidence 
of the occurrence of drag reduction in absence of boundary layers. A weakly non-linear model 
proposed by Fermi for the growth of the mixing layer is reported in the Appendix. 



1 



I. INTRODUCTION 



Polymer additives have dramatic effects on the dynamics of turbulent flows, the most 
important being the reduction of turbulent drag up to 80 % when few parts per million of 



long-chain polymers are added to water The paramount relevance of this phenomenon 
motivated the strong efforts of researchers aimed to achieve a better understanding of the 
basic mechanisms of polymer drag reduction. The natural framework of drag-reduction 
studies is the case of pipe flow, or channel flow. Within this context the reduction of the 
frictional drag against material wall originated by the addition of polymers, manifests as an 
increase of the mean flow across the pipe or channel at given pressure drop. 

Recent studies (see, e.g., showed that a drag-reduction phenomenon may also 

occurs in the absence of physical walls. In this case the drag which is reduced is not the 
frictional drag against the boundaries of the flow, but the turbulent drag of the bulk flow 
itself. In particular, in the case of homogeneous isotropic turbulence, it has been observed a 
reduction of the rate of energy dissipation at fixed kinetic energy associated with a reduction 
of velocity fluctuations at small scales^- 6|. In turbulent systems with a non- vanishing local 
mean flow (e.g the Kolmogorov flow), it has been shown that polymers causes a reduction of 
the Reynolds stresses which results in an increased intensity of the mean velocity profile j2|. 
This phenomenon, which occurs in absence of boundaries, is remarkably similar to increase of 
throughput observed in pipe or channel flows, and suggest the existence of common features 
and possibly of common physical mechanisms between the drag-reduction occurring in wall- 
bounded and in bulk flows. 

In the present paper we provide further evidence of turbulent drag reduction in bulk 
flows by studying the effects of polymers additives in the Rayleigh- Taylor (RT) setup of 
turbulence convection. A previous study [7] has already shown that polymers affect the 
early stage (linear phase) of the RT instability which occurs at the interface between two 
unstably stratified fluids. Here, we show that polymers also induce strong modifications in 
the dynamics of the turbulent mixing layer, which develops in the late stage of the mixing 
process. In particular we study how polymers are able to affect the process of turbulent heat 
transfer with a mechanism which is probably more general than the particular case studied 
here. 

Preliminary results have been presented in js] and are briefly reported here for com- 



pleteness. We provide here new results supporting our interpretation of the mechanism at 
the basis of the observed effects together with results on polymer statistics and small scale 
turbulence. 

The remaining of this paper is organized as follow. In Section |TT] we introduce the vis- 
coelastic Rayleigh-Taylor problem and give some details on the numerical strategy we ex- 
ploited to study polymer dynamics. In Section IIIII we analyze the statistics of polymer 
elongations. In Section IIVI we show the effects of polymers on the turbulent mixing. In 
Section |V] we discuss the drag reduction phenomenon in the viscoelastic RT. In Section I VI I 
we study the effects induced by polymers on the heat transport. Conclusions are devoted 
to a short discussion on the possibility to observe the described effects in the laboratory. 
Finally, in the Appendix we briefly describe the model for the growth of the mixing layer 
proposed by Fermi. 



II. THE VISCOELASTIC RAYLEIGH TAYLOR MODEL 

We will focus our attention on the miscible case of the RT system at low Atwood number 
and Prandtl number one. Within the Boussinesq approximation, generalized to a viscoelastic 
fluid using the standard Oldroyd-B modeljo], the equations for the dynamics of the velocity 
field u coupled to the temperature field T(x, t) (which is proportional to the density via the 
thermal expansion coefficient (3 as p = po[l — f3(T — Tq)], po and Tq are reference values) and 
the positive symmetric conformation tensor of polymers aij{:s.,t) read: 

+ u • Vu = -Vp + z/V^u - (3gT + —V ■ cr 



dtT + u-VT = kV^T (1) 
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dta + w-Vo = (Vu)^ • a + (7 ■ (Vu) (a - I) + fi;„VV 

together with the incompressibility condition V ■ u = 0. In ([T]) g = (0, 0, —g) is gravity 
acceleration, v is the kinematic viscosity, n is the thermal diffusivity, rj is the zero-shear 
polymer contribution to total viscosity vt = + 1]) (proportional to polymers concentra- 
tion) and Tp is the (longest) polymer relaxation time, i.e. the Zimm relaxation time for a 
inear chain Tp = vR^/ {pksT) with ks Boltzmann constant and Rq the radius of gyration 
^. The diffusive term KpV^cr is added to prevent numerical instabilities jlol |. 

The initial condition for the RT problem is an unstable temperature jump T(x, 0) = 
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512 1024 27r An 1 0.5 3 • 10 
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512 1024 2tt An 1 0.5 3 • lO"'' lO^^ 0.2 1 



B 



512 1024 2tt An 1 0.5 3 • lO"'' lO^^ 0.2 2 



C 



512 1024 2n An 1 0.5 3 • lO"'' lO^^ 0.2 10 



TABLE I. Parameters of the simulations 



— {6o/2)sgn{z) in a fluid at rest u(x, 0) = with coiled polymers cr(x, 0) = I. The physi- 
cal assumptions under which the set of equations ([T]) is valid are of small Atwood number 
A = {1/2)P6q and dilute polymers solution. Experimentally, density fluctuations can also 
be obtained by some additives (e.g., salt) instead of temperature fluctuations: within the 
validity of Boussinesq approximation, these situations are described by the same set of equa- 
tions ([I]). In the following, all physical quantities are made dimensionless using the vertical 
side, Lz, of the computational domain, the temperature jump 9o and the characteristic time 
r = (Lz/ AgY^'^ as fundamental units. 

Numerical simulations of equations (II]) have been performed with a parallel pseudospec- 
tral code, with 2-nd order Runge-Kutta time scheme on a discretized domain of N^xNyX 
grid points. Periodic boundary conditions in all directions are imposed. The initial pertur- 
bation is seeded in both cases by adding a 10% of white noise (same realization for both 
runs) to the initial temperature profile in a small layer around the instable interface at 2; = 0. 
Because of periodicity along the vertical direction, the initial temperature profile has two 
temperature jumps: an unstable interface at 2; = which develops in the turbulent mixing 
layer and a stable interface at at 2; = -^Lz/2. Numerical simulations are halted when the 
mixing layer is still far from the stable interface, whose presence has no detectable influence 
on the simulations (velocities there remain close to zero) . The results of the reference New- 
tonian simulation (denoted by run A^) are compared with those of three viscoelastic runs 
[A^B and C) with identical parameters and different polymer relaxation time (see Table [Tj). 

III. STATISTICS OF POLYMER ELONGATIONS 

Before presenting the results of our numerics, let us discuss the theoretical behavior 
expected for polymers statistics in the "passive case" in which their feedback on the flow 
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is neglected. Recent studies of Newtonian RT turbulence |llHl3| support the picture of a 
Kolmogorov scenario, in which the buoyancy forces sustain the large scale motion, but they 
are overcome at small scales by the turbulent cascade process. The accelerated nature of the 
system results in an adiabatic growth of the flux of kinetic energy in the turbulent cascade 
e ~ {AgYt. As a consequence, the Kolmogorov viscous scale i] ~ i/'^/'^e~^^'^ and its associated 
time-scale ^ (zz/e)^/^) decrease in time as 77 ~ iy^^^(Ag)~^^'^t~^^'^ and ~ v^^'^{Ag)^^ ~ 
^-1/2 respectively. 

The Weissenberg number Wi = Tp/r^, which measure the relative strength of stretching 
due to velocity gradients and polymer relaxation, grows as Wi ~ Therefore, even if the 
relaxation time of polymer Tp is sufficiently small to keep the polymers in the coiled state in 
the initial stage of the evolution, they are expected to undergo a coil-stretch transition as 
the system evolves. The Lumley scale, defined as the scale ii whose characteristic time is 
equal to the polymer relaxation time r^^ ~ — Tp grows in time as £l ~ AgTp^'^t^^'^. 

In view of the fact that the turbulent inertial range extends from the integral scale C ~ Agt'^ 
to the dissipative scale i] ~ t^^^^ the temporal evolution of the Lumley scale guarantees that 
for long times one has r] < < C 

It is worth noting that in two dimension the behavior would be the opposite. In con- 
trast to the three-dimensional (3D) case, the phenomenology of RT turbulence in 2D is 
characterized by a Bolgiano scenario, which originate from a scale-by-scale balance between 
buoyancy and inertial forces . The resulting scaling behavior of velocity increments 

is 6iU ^ {Agy^^^i^^^t~^^^ , which gives for the dissipative scale rj ~ {Ag)^^^'^i'^^^t^^^ and 

~ (Ag)^^^'^i'^^^t^^^. Therefore in the 2D case the Weissenberg number decreases in time 
as Wi ~ and polymers will eventually recover the coiled state. The Lumley scale 

decay as £l — AgTp^'^t~^/'^ , and in the late stage of the evolution will become smaller that 
the dissipative scale 77. 

Under the hypothesis that these scaling behaviors remain valid also in the presence of 
polymer feedback to the flow, one may conjecture that viscoelastic effects in 3D RT turbu- 
lence become more and more important as the system evolves (while, as explained, in the 2D 
case they are expected to be transient, and to disappear at the late stage of the evolution). 

The presence of a coil-stretch transition in the 3D RT flow is confirmed by the behavior 
of the rms polymer elongation Rrms = (('^^(<^)/3))^^^-Ro measured in our simulations (see 
insets of Figs. [H El [3]). In the initial stage of the evolution the velocity gradients are too 
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FIG. 1. Pdfs of polymers elongation at times 3 < t/r < 4. Solid line t/r = 3.1, dashed line 
t/r = 3.4, dotted line t/r = 3.7, dash-dotted line t/r = 4.0. Inset: Rms of polymers elongation as 
a function of time. Data from run A. 




FIG. 2. Pdfs of polymers elongation at times 2.5 < t/r < 3.5. Solid line t/r = 2.5, dashed line 
t/r = 3.8, dotted line t/r = 3.1, dash-dotted line t/r = 3.4. Inset: Rms of polymers elongation as 
a function of time. Data from run B. 

weak to significantly stretch the polymers, and Rrms ~ -Ro- At time t c:^ t a transition 
occurs, and polymers start to elongate. After a transient exponential growth, a regime 
characterized by a linear growth Rrms ~ t sets in, which is consistent with the growth of 
elastic energy discussed in Section |Vl The pdf of elongations in this stage of the evolution 
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FIG. 3. Pdfs of polymers elongation at times 2.2 < t/r < 3.2. Solid line t/r = 2.3, dashed line 
t/r = 2.5, dotted line t/r = 2.8, dash-dotted line t/r = 3.1. Inset: Rms of polymers elongation as 
a function of time. Data from run C. 

are not stationary, but their right tail collapse once rescaled with Rrms (see Figs. [H Ej |3]). 
Oldroyd-B model allows a priori for infinite elongations, but we observe an exponential cutoff 
for the right tail of the pdfs, which is a genuine viscoelastic effect: polymer feedback is able 
to reduce the stretching efficiency of the flow. These observations lead to the conclusion that 
polymers dynamics follows adiabatically the accelerated growth of the flow and generates a 
strong feedback which manifests in the appearance of a cutoff for their elongations. 




IV. EFFECTS OF POLYMERS ON MIXING PROPERTIES 



The evolution of the turbulent mixing layer is strongly affected by polymer additives. 
For a Newtonian flow, because of the constant acceleration provided by the gravity force, 
one expects the width h{t) of the mixing layer to grow as h{t) = aAgt'^, where a is a 



dimensionless parameter to be determined empirically 15Nl7l|. Several definitions of h{t) 



have been proposed, based on either local or global properties of the mean temperature 



profile T{z,t) (the overbar indicates average over the horizontal directions) [12|, |18N20|. 
Here, we adopt the simplest measure hr based on the threshold value of z at which T{z,t) 
reaches a fraction r of the maximum value i.e. T{±hr{t)/2,t) = ^r6Q/2. 



In the viscoelastic solution the growth of the mixing layer is faster than in the Newtonian 
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FIG. 4. Temporal evolution of the mixing layer width h{t). Inset: fraction of mixed fluid within 
the mixing layer p = \/{LxLyLz) f d'^x9[{r9Q/2)'^ — T^]. Newtonian flow: solid line. Viscoelastic 
flows: dashed line (Run A), dotted line (Run B), dash-dotted line (Run C). 
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FIG. 5. Time evolution of the ratio Wrms/urms (panel a) and velocity gradients (panel b). Correla- 
tion length of vertical {Rw, panel c) and horizontal {Ru, panel d) velocity components. Newtonian 
flow: solid line. Viscoelastic flows: dashed line (Run A), dotted line (Run B). 

case (see Fig. H]), and the acceleration effect is stronger for polymers with longer relaxation 
times. On a coarse scale this effects produces a mixing enhancement. On the other hand, 
the viscoelastic fluid is less uniformly mixed within the mixing layer itself. In the Newtonian 
case the volume of the region where |T(x, t)| < r9o/2 is roughly 80% of the volume of the 
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mixing layer at the same time. Conversely, the fraction of mixed fluid within the mixing 
layer reduces up to 50% for the viscoelastic case (see inset of Fig. H]). These results indicate 
that the effects of polymers on the mixing efficiency is twofold. At large scale they enhance 
the mixing by accelerating the growth of the mixing layer and at small scale they reduce 
the mixing efficiency of the turbulent flow. 

These effects are accompanied by an increase of the anisotropy of the flow. In Figure [5] we 
show the ratio between rms of vertical {wrTns) and horizontal velocities (u^ms) and velocity 



gradients. The velocity ratio, which is around 1.8 for the Newtonian case [13|], becomes larger 
than 2.5 for the viscoelastic run. This phenomenon is associated with the enhancement of 
the vertical velocity of the mixing layer. The reduction of small-scale mixing efficiency 
results in the persistence of anisotropy also at small scales (i.e. in the velocity gradients), 
at variance with the Newtonian case in which it is almost absent. The viscoelastic flow is 
therefore characterized by the presence of faster and larger plumes than those characterizing 
the Newtonian case. The increased coherence of thermal plumes can be quantified in terms 



of the enhancement of the ve 



ocity correlation length (here defined as the half width of the 



velocity correlation function 12|), of both horizontal and vertical velocity components (see 
Fig.El). 

V. INTERPRETATION IN TERMS OF DRAG REDUCTION 

The energy balance of the viscoelastic RT system differs from the Newtonian case because 
of the elastic contribution to the energy and dissipation. The energy can be written as the 
sum of potential, kinetic and elastic contributions: 

E = P + K + J: = -l3g{zT) + hu^) + ^ [{tra) - 3] (2) 

2 Tp 

and the energy balance reads: 

dE 

^ = (3) 

where Si, = ^{{daUpY) is the viscous dissipation and the last term represents elastic dissipa- 
tion Ey. = TLjTp. The evolution of the system is sustained by the consumption of potential 
energy, which provides a power source —^ = f3g{wT) (where w is the vertical velocity 
component). It is worth noting that the rate of energy injection is not determined a priori. 
Indeed, it is the dynamics itself which determines the rate of conversion of potential energy 
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FIG. 6. Temporal evolution of the potential energy P (panel a), kinetic energy K panel (b), 
elastic energy S (panel c) and viscous energy dissipation (panel d). Newtonian flow: solid line. 
Viscoelastic flows: dashed line (Run A), dotted line (Run B). 



into kinetic and elastic energy. Our numerics reveals that polymers accelerate this process 
(see Fig.E]), and that kinetic energy for viscoelastic runs is larger than that of the Newtonian 
case (of about 40% at t = 3.5r). We remark that the faster growth of kinetic energy is not 
a straightforward consequence of the speed-up of potential energy consumption, due to the 
accelerated growth of mixing layer. Part of the potential energy is indeed converted into 
elastic energy and finally dissipated by polymers relaxation to equilibrium. 

The increase of kinetic energy is accompanied by a reduction of viscous dissipation 
(Fig. Ell). This is a clear fingerprint of a drag reduction phenomenon as defined for 
homogeneous-isotropic turbulence 0, 0], i.e. a reduction of turbulent energy dissipation 
at given kinetic energy. In the present case, a quantitative measure of the drag reduction 
is provided by the ratio between the loss of potential energy and the resulting plumes ki- 
netic energy. The first can be easily computed by the definition of the potential energy 
P = —Pg{zT) assuming a linear temperature profile within the mixing layer, which gives 
AP = P{0)—P{t) ~ l/6Agh{t). An estimate of the kinetic energy associated with large scale 
plumes can be obtained in terms of the mixing layer growth rate h{t) as ~ l/2(/i(t))^. 
We remark that a similar estimation was proposed by Fermi for modeling the growth of 
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FIG. 7. Time evolution of the drag reduction factor /. Newtonian flow: solid line. Viscoelastic 
flows: dashed line (Run A), dotted line (Run B), dash-dotted line (Run C). 
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FIG. 8. Energy spectra of the vertical velocity component at time t = 3.1r Newtonian flow: solid 
line. Viscoelastic flows: dashed line (Run A), dotted line (Run B), dash-dotted line (Run C). 



mixing layer (see Appendix). The drag reduction coefficient / is then defined as 



(4) 



Kl ' Via 

which turns out to be inversely proportional to the coefficient a which characterizes the 
mixing layer growth rate 20||. With this definition, we measure 22% of drag reduction for 
the viscoelastic run B and 30% for the run C (see Fig. [7j) 
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The scenario which emerges from these results is that polymers reduce the turbulent drag 
between rising and sinking plumes. The RT viscoelastic system is therefore able to convert 
more efficiently potential energy into kinetic energy contained in large plumes. Conversely, 
the turbulent transfer of kinetic energy to small-scale structures is reduced, which results in 
a reduction of the viscous energy dissipation. This picture is confirmed by the inspection of 
the energy spectra (see Fig. [8]). At small scales we found a suppression of turbulent kinetic 
energy with respect to the Newtonian case, while at large scale an increase of the kinetic 
energy is observed. 

The drag reduction scenario therefore provides a clear interpretation of the effects ob- 
served on the mixing properties. The enhancement of large-scale mixing associated with 
the faster growth of the mixing layer is directly connected to the reduced friction between 
plumes, and the reduced efficiency of small-scale mixing is a natural consequence of the 
suppression of small-scale turbulence. 

The accelerated nature of the RT turbulence poses an interesting question about the exis- 
tence of an asymptotic state for viscoelastic RT. For the Newtonian case the phenomenolog- 
ical theory assumes that in the late stage of the evolution all terms in the energy balance ([3]) 
have the same temporal scaling determined by gravity forces. This implies that ^ ^ t 
and K t^. In the viscoelastic case it is not possible to fix a priori the scaling law for the 
elastic contribution, because elastic energy S is proportional to the elastic dissipation rate 
e-z = 2S/rp. Assuming that the latter has the same temporal scaling than the viscous dissi- 
pation, ~ ~ t, one gets that the elastic contribution to the total energy should become 
negligible at long times. On the other hand, the assumption that elastic and kinetic energy 
have the same scaling E ~ K ~ leads to the conclusion that elastic dissipation would 
eventually dominate over the viscous one. Our simulations support the second hypothesis: 
the ratio between elastic end viscous dissipation is not constant, and grows almost linearly 
in time (see Fig. [9]). A deeper investigation of this asymptotic state in which polymers are 
strongly elongated would require to go beyond Oldroyd-B model, and to adopt more realistic 
polymer model (e.g. FENE-P model) which accounts for maximal elongation and non-linear 
relaxation. 
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FIG. 9. Temporal evolution of the ratio between elastic and viscous dissipation. Run A: dashed 
line. Run B: dotted line. 
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FIG. 10. Temporal evolution of Nusselt number Nu = {wT)h/ {ndo). Inset: Nusselt number vs. 
Rayleigh number Ra = Agh? /{uk). Newtonian flow: solid line. Viscoelastic flows: dashed line 
(Run A), dotted line (Run B), dash-dotted line (Run C). 

VI. HEAT TRANSPORT ENHANCEMENT 

The heat transport efficiency in turbulent convection is usually measured by the Nusselt 
number Nu — {wT)h/ {k6q), which represents the ratio between convective and conductive 
heat transport. For a developed turbulent flow the Nusselt number is expected to behave as 
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FIG. 11. Temporal evolution of the contributions to the heat transfer efficiency. Panel (a): Mixing 
layer width h{t). Panel (b): Rms vertical velocity Wrms- Panel (c): Rms temperature T^ms- Panel 
(d): Correlation between temperature and vertical velocity CwT{z^t). Newtonian flow: solid line. 
Viscoelastic flows: dashed line (Run A), dotted line (Run B). 



a simple power law with respect to^the dimensionless temperature jump which defines the 
Rayleigh number Ra = Agh^ /{vk) [21]. For a flow in which boundary layers are irrelevant, 
as in our case, Kraichnan |22| predicted many years ago an asymptotic regime which is 
expected to emerge at very large Ra. For this so-called ultimate state regime of thermal 
convection dimensional analysis predicts the scaling laws 21] 



Nu ~ Pr^/^Ra^/^ Re ~ Pr-^l'^Ra}''^ . 



(5) 



For the case of time dependent RT turbulent convection, all these dimensionless quan- 
tities depend on time. Dimensionally estimation gives (for the Newtonian case) Ra ~ 
{(3g9Q)H^/{h'K), Re ^ {f3g9o)H^/h' and Nu ~ {(3g9o)H^/K, which indeed imply the scaling 
aws (|5]) and which have been observed recently in numerical simulation of RT turbulence 
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The addition of polymers strongly enhances the efficiency of heat transport, i.e. the 
Nusselt number grows faster both as a function of time and as a function of Ra jsj, and the 
effects increase with the polymer relaxation time, as shown in Fig. [TOl In order to identify 
the different causes which contribute to this effect, it is useful to rewrite the Nusselt number 
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as: 

Nu = -\-hWrmsTrmsCu,T (6) 

where Cu)T = (wT) / (wrmsTrms) is the correlation between the vertical velocity component 
and the temperature field. In the four panels of Fig. [TT] we plot the four contributions h 
(panel a) Wrms (panel b), Trms (panel c) and C^t (panel d) as a function of time. It is 
evident that the increased heat transfer is not simply a consequence of the faster evolution 
of the mixing layer h, but also of the increased rms of the vertical velocity component, 
Wrms- Moreover, the reduction of small-scale turbulent mixing causes an increase of the 
temperature fiuctuations Trms which also gives a positive contribution to the Nusselt number. 
Finally, in the viscoelastic case we found stronger correlations between temperature and 
vertical velocity component which therefore transport heat more efficiently. In conclusion, 
the increased heat transport efficiency is a combined effect of different contributions: the 
presence of faster thermal plumes, the reduced turbulent mixing, and the stronger correlation 
between thermal plumes and the vertical velocity component. While the first contribution 
is distinctive of RT turbulence, the others could in principle be observed in other thermal 
convective systems. A recent experiment performed within the framework of Rayleigh- 
Benard convection indicates in that case an opposite effect of heat transfer reduction 24 1 
but this can be probably attributed to the moderate stretching of polymers in that case. 



VII. CONCLUSIONS 



The behavior of viscoelastic flows in the RT setup provides the flrst clear evidence of 
simultaneous occurrence of both polymer drag reduction and heat transport enhancement. 
Drag reduction in this system is caused by a reduced drag between rising and sinking thermal 
plumes, a fact which implies the speed up of the mixing layer growth. This process shares 
many analogies with drag reduction observed in homogeneous, isotropic turbulence, namely 
the suppression of small-scale turbulence which results in a reduced viscous drag. These 
analogies provide a support to the conjecture of a common underlying mechanisms behind 
these different manifestations of the polymer drag reduction in bulk flows. 

For RT system it is possible to introduce a drag coefficient in terms of the ratio between 
the potential energy loss which forces the flow, and the resulting kinetic energy associated 
with thermal plumes. The viscoelastic case is characterized by faster and more coherent 
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thermal plumes. The effects on mixing is to enhance the large-scale mixing, and to reduce 
the small-scale one. As a consequence, the drag coefficient is reduced and the heat transport 
efficiency, measured by the Nusselt number, is increased. 

We conclude with some speculations on the possible observation of heat transfer enhance- 
ment in laboratory experiments. The values of the parameters used in our simulations can be 
used to determine the setup for a comparable experiments. The units of time T and length 
L which allow one to convert the parameters of our simulations into physical quantities can 
be fixed by matching the numerical values of viscosity D — ?> ■ 10~^ and gravity g — (4A)~^ 
used in our simulations with physical values g — 9.81ms~^, v — vh^o — 10""^m^s~^: 

,2 



I>2 g 



(7) 



T'^t-, (8) 
V g^ 

By choosing the Atwood number A = 0.1 one gets L ~ 1.4cm and T ~ 0.06s. This 
correspond to an experimental box of Lx,y — 10cm and — 20cm, and polymer relaxation 
times Tp = 60ms for the case A, which is close to realistic relaxation times of long-chain 
polymers in water. The evolution of the system will be quite fast: the time required for 
the mixing layer to invade the whole box is estimated to be roughly 2s. Let us notice that 
the limit of small Atwood number, required in the present study to justify the Boussinesq 
approximation, is not a constraint for an experimental setup, where large values of A can 
be obtained by means of some additives (e.g. salt) to generate density differences. It would 
be interesting to observe experimentally the influence of non-Boussinesq effects on the drag 
reduction phenomenon. 
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Appendix A: Fermi model for the growth of mixing layer 



Enrico Fermi, together with John von Neumann, were probably the first who considered 
a model for the growth of mixing layer in the nonlinear stage. The model is described in 
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FIG. 12. Fermi's model for the evolution of the interface. The smooth interface ACE' is replaced 
by the square wave ABCEE' . 

two reports of the Los Alamos Scientific Laboratory, the first from September 1951 (Fermi 
alone) and the second from August 1953 (Fermi and von Neumann) never published 25 |. 
The idea of this work, in the words of the authors, is to "discuss in a very simplified form 
the problem of the growth of an initial ripple on the surface of an incompressible liquid in 
presence of an acceleration". The first report of Fermi considers the interface between a 
liquid and vacuum, while the second report with von Neumann analyzes the case of two 
fluids of different densities. 

The idea of Fermi is to approximate the interface with a square wave whose shape is 
characterized by three parameters: the heights of spike and bubble a and h and the width 
of the spike x (see Fig. [T2|) . Incompressibility gives a relation among these quantities, 
h = ax/{l — x). Fermi next considers the Euler-Lagrange equations for the variation of the 
potential and kinetic energy and obtains a couple of equations for the evolution of a and x. 
In the following we consider a simplified version of Fermi model with bubble-spike symmetry 
{h = a, X = 1/2), consistent with the Boussinesq approximation discussed in the present 
paper. 

The variation of potential energy to generate the profile in Fig. [12] is 



U 



— 7^ gLxLya^ 



(Al) 



For the kinetic energy, assuming that the "plumes" ABCO and CC'E'E move respectively 
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up and down with velocity a and plumes BB'C'C and OCED move respectively right and 
left with the same velocity (for incompressibility) one obtains 

K = Pl^J^L^Lyac? (A2) 

From the Lagrange equations 

_ dl^ _ _dU_ 
dt da da da 

one the obtains a differential equation for the a{t) without free parameter. We remark that 
[3]) assumes that all potential energy is transformed in large scale kinetic energy generating 



the motion of the interface. In a later stage, in which a turbulent flow develops, we can still 
try to use (]A3ll but now with a factor < 5 < 1 in front of the rhs which takes into account 
that a fraction (1 — 5) of potential energy goes in viscous dissipation (through the turbulent 
cascade). With this correction the equation for a{t) reads 

cPa 1 / da\^ 1 ^ ^ . ^ . 

Introducing the width of the interface h{t) = 2a{t) and solving (1A4I) with initial condition 
h{0) = ho and replacing 5 = 80; we finally get 

h{t) = ho + {AaAghoY^^t + aAgf (A5) 

which is the form proposed from other authors on the basis of completely different consid- 
erations |20|, 126 1 . 
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